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SUMMRY 

^t  is  shown  that,  for  ideal  bundles  having  a  given  number  of  fibers  that 
will  break  under  a  load  per  fiber  Zp  and  a  given  number  of  fibers  that  will 
not  break  under  a  load  Zg  (^^l)^  "tbe  average  number  of  unbroken  fibers  in  the 
bundle  at  loads  between  Z^  and  Zg  depends  linearly  on  the  probability  that 
an  individual  fiber  picked  at  random  will  support  a  load  Z.  Bundles  are 
grouped  according  to  the  number  of  fibers  in  the  bundle  that  can  support  that 
load  per  fiber  at  which  the  average  total  load  supported  by  the  bundle  is  a 
maximum.  An  approximation  to  the  maximum  of  the  average  load  supported  by  the 
bundle  is  arrived  at  for  each  group,  and  this  maximum  averaged  over  the  groups 
constitutes  a  lower  bound  for  the  average  breaking  strength  of  the  bundle.  The 
results  agree  with  those  of  Pierce  and  Daniels^^ 


INTRODUCTION 

A  classical  fiber  is  one  that  will  support  a  load  less  than  its  breaking 
strength  indefinitely  without  stretching  or  breaking  but  will  break  immediately 
under  any  load  equal  to  or  in  excess  of  its  breaking  strength.  The  fundamental 
papers  dealing  with  the  breaking  strengths  of  an  ideal  bundle  (no  friction  or 
twist)  of  equal-length  classical  fibers  are  those  by  Pierce  (ref.  l)  and  Daniels 
(ref.  2).  Pierce  was  the  first  to  obtain  an  approximation  to  the  average  break¬ 
ing  strength  of  a  bundle  of  N  such  fibers.  Daniels  rigorously  showed  by 
sophisticated  mathematical  methods  that  the  expression  of  Pierce  is  the  correct 
asymptotic  limit  of  the  breaking  strength  as  N  gets  large.  He  also  found  the 
asymptotic  limit  of  the  variance  of  the  breaking  strength  for  large  N.  In 
this  report,  the  method  of  reasoning  employed  leads  to  the  establishment  of  the 
aforementioned  results  from  a  viewpoint  that  is  based  on  physical  rather  than 
mathematical  insight. 


SYMBOLS 

B(Z)  probability  that  randomly  selected  fiber  will  break  at  load  less  than  Z 

b(Z)  frequency  function  for  breaking  strength  of  individual  fibers  randomly 
selected 

E  expectation  value 


G(  z )  mil) 


G]^  derivative  of  G(  Z )  with  respect  to  Z  evaluated  at  Z 

Hn(y)  Hermite  polynomial  of  y 

k  dujnmy  index 

L  total  load  on  bundle  of  fibers 

L( I )  total  load  at  I  averaged  over  population  of  bundles 
I  load  per  fiber 

Z  load  per  fiber  at  which  average  total  load  L  is  maximum 

N  number  of  fibers  in  bundle 

n  number  of  unbroken  fibers  in  bundle 

n  number  of  unbroken  fibers  in  bundle  at  load  Z 

0  order 

P  probability 

R(Z)  probability  that  randomly  selected  fiber  will  support  load  up  to  Z 

5  average  breaking  strength  of  bundle 

y  variable  (li[R(Z)  -  n) 

&£  (mil)  -  n)/(N  -  n) 

6  variable  of  integration 

A£  value  of  I  that  maximizes^average  load  supported  by  bundles  having  n 

unbroken  fibers  at  load  I 

third  moment  about  origin  for  variable  y 

a  standard  deviation 

CONCEPT  OF  PATH 

Let  the  distribution  of  breaking  loads  of  a  classical  fiber  be  given  by 
the  frequency  function  b(Z)^  where  b(Z)di  is  the  probability  that  a  randomly 
selected  individual  fiber  will  have  a  breaking  strength  lying  in  dZ  at  Z. 
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The  cumulative  distribution  function 


B(Z)  =  f  b(0)d0 

is  the  probability  that  a  randomly  selected  fiber  will  break  at  some  load  less 
than  Z .  Let  R( Z )  =  1  -  B( Z )  denote  the  probability  that  a  fiber  will  remain 
unbroken  up  to  load  I . 

When  a  load  L  is  hung  on  an  ideal  bundle  of  W  classical  fibers^  the 
load  distributes  itself  evenly  over  the  unbroken  fibers  so  that  each  of  the  n 
unbroken  fibers  supports  a  load  I  =  L/n.  The  breaking  strength  of  a  given 
fiber  of  the  bundle  does  not  depend  on  the  action  of  the  other  fibers  of  the 
bundle  and^  hence ^  the  breaking  strength  of  the  bundle  may  be  determined  from 
the  strengths  of  the  individual  fibers.  In  particular,  suppose  one  considers 
the  bundle  as  separated  into  its  independent  individual  fibers  with  separate 
loads  I  of  equal  magnitude  on  each  fiber.  Let  this  load  I  be  allowed  to 
increase  monotonically  from  zero  to  the  breaking  strength  of  the  strongest 
fiber.  Then,  at  any  point  in  the  process  where  the  individual  fibers  support  a 
load  Z,  the  total  load  L(Z)  supported  by  the  bundle  is  the  product  of  the 
number  of  unbroken  fibers  n( Z )  and  the  load  Z  acting  on  each  unbroken  fiber . 
Thus,  the  breaking  strength  of  the  bundle  is  the  maximum  of  the  total  load  L. 
After  this  maximum  L  has  been  reached,  further  increases  in  Z  cause  a  reduc¬ 
tion  in  the  total  load  L(Z)  =  Zn(Z).  (if  the  fibers  were  not  separated  in  this 
manner,  and  if  breaking  strength  were  determined  by  increasing  L  until  break¬ 
age  occurred,  no  decreases  in  total  load  L  would  appear,  and  the  bundle  would 
immediately  break  after  the  breaking  strength  of  the  bundle  had  been  reached.) 

The  situation  may  be  depicted  by  means  of  figure  1  (p.  4).  In  fig¬ 
ure  1(a),  a  typical  plot  of  Z  against  R( Z)  is  given.  (Note  that  in  the  fig¬ 
ures,  R(Z)  decreases  from  unity  to  zero  as  one  proceeds  from  left  to  right  along 
the  horizontal  axis. )  It  can  be  seen  that  as  Z  increases  from  zero  to  high 
values,  R(z)  decreases  from  unity  to  zero  or,  if  R(z)  is  to  be  considered  the 
independent  variable,  Z  increases  from  zero  to  high  values  as  R(Z)  decreases 
from  unity  to  zero.  The  path  followed  by  a  bundle  of  N  fibers  averaged  over 
the  population  of  bimdles  is  shown  as  the  straight  line  in  figure  1(b),  where 
the  expected  number  of  surviving  fibers  is  given  by  E(n)  =  NR(Z).  The  path 
taken  by  a  typical  bundle  is  shown  in  figure  1(c),  where  n  for  that  bundle 
is  plotted  against  R(z).  For  large  N  at  least,  the  actual  path  followed  by 
a  particular  bundle  would  be  expected  to  lie  fairly  close  to  the  average  path 
as  given  by  the  E(n)  line. 


RESULTS  OF  PIERCE  AND  DANIELS 

The  average  load  supported  by  the  bundle  as  a  function  of  Z  is 

L(Z)  =  ZE(n)  =  NZR(Z)  (l) 

The  maximum  average  load  supported  by  the  bundle  is 

=  L(Z)  =  NZR(Z)  (2) 
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(a)  Load. 


(b)  Expectation  of  survival. 


Probability  that  fiber  will  support  load  I,  R(0 
(c)  Number  of  unbroken  fibers. 

Figure  1.  -  Plots  with  R{1)  as  independent  variable. 


where  I  satisfies  the  equation 


dL(Z) 

dl 


R(Z)  =  Zb(z) 


Equations  (2)  and  (s)  represent  the 
approxiiriation  to  the  breaking  strength 
of  the  bundle  as  given  by  Pierce. 
Daniels  showed  that  this  result  is  the 
correct  asymptotic  limit  for  the  ex¬ 
pectation  value  of  the  breaking 
strength  for  large  N.  He  also  found 
that  the  asymptotic  limit  for  the 
standard  deviation  of  the  strength  for 
large  N  is  given  by 


0  =  1 


Daniels  makes  the  pertinent  remark  that 
equation  (4)  would  follow  if  one  could 
assume  that  the  breaking  strength  is^ 
only  dependent  on  the  number  n^=  n( Z ) 
of  fibers  surviving  at  a  load  Z  where 
n  is,  of  course,  distributed  according 
to'  the  simple  binomial  law.  He  adds, 
however,  that  there  appears  to  be  no 
”a  priori  justification"  for  this 
assumption. 

The  average  breaking  strength  S 
of  a  bundle  of  N  fibers  is  the  maxi¬ 
mum  load  that  the  bundle  supports 
averaged  over  the  population  of  bundles. 
For  a  particular  bundle,  this  maximum 
load^will  not,  in  general,  occur  at 
Z  =  Z,  and  hence,  it  will  be  larger 
than  L( Z ) ,  the  value  of  L  ^at  Z .  On 
the  other  hand  I^ax  L( Z )  averaged 

over  the  population  of  bundles,  and 
therefore,  Lj^ax  constitutes  a  lower 


_  therefore,  Lj^ax  constitutes  a  lower 

bound  for  S.  If  the  maximum  load  for  every^possible  path  (in  the_sense  of 
fig.  1(c))  occurred  at  Z  =  Z,  then  Lmax  =  L( Z )  would  indeed  be  S  exactly. 
It  can  be  said  that  the  approximation  to  S  as  given  by  equation  (2)  treats 
every  path  as  though  its  maximum  load  occurred  at  z  =  ?. 


AVERAGE  PATH  BETWEEN  TWO  DESIGNATED  POINTS 


It  shall  first  be  shown  that,  for  values  of  Z  lying  between  Zq  and 
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the  average  path  taken  by  bundles  which  are  constrained  to  go  through  the 
points  R(Zx)^  n^  and  R(Z2)^n2  (i.e.^  bundles  which  have  unbroken 
fibers  at  a  load  l\  and  ns  unbroken  fibers  at  a  load  Zg)  is  merely  the 
straight  line  connecting  these  two  points  on  a  plot  of  n  against  R( Z ) .  The 
proof  is  straightforward  and  goes  as  follows. 

The  probability  that  a  bundle  originally  consisting  of  N  fibers  takes 
a  path  such  that  it  has  np  survivors  at  Zp  and  ns  survivors  at  Zs 
where  ^  >  ng  >  0  and  0  <  <  ^S 

P(n']  at  Zn^  no  at  Zo)  =  ^ - .......El - - - 

^  ^  (1^  -  np)I(ni  -  ii^nsl 


[B(z^)]N-ni[R(j^)  .  R(i2)]ni-n2[R(^2)]^2  (5) 

The  probability  that  a  bundle  originally  consisting  of  N  fibers  takes  a  path 
that  passes  through  the  three  points  np  at  Zp^  n  at  Z^  and  ns  at  Zs^ 
where  ^>3a2_>n>n2>0  and  0  <  ^p  £  ^  <  ^S^ 

P(ni  at  li,  n  at  I,  nz  at  Iz)  . 

[b(Ii)]''-1[k(*i)  -  [e(1)  -  (6) 

The  conditional  probability  that  a  bundle  which  passes  through  np  at  Zp  and 
n2  at  Z2  has  n  unbroken  fibers  at  Z  is  given  by  dividing  equation  (6) 
by  equation  (s) 

( np  -  n2 ) I 

P(n  at  Z/ni  at  ng  at  Zg)  =  (  — V'nyi^n  - 

[r(Zi)  -  R(Z)]^^~^[r(Z)  -  R(Z2)]^'^^ 

[r(Zi)  -  R(Z2)]"1-"2 

The  expectation  value  of  n  at  Z  for  a  bundle  passing  through  np  at  Zp 
and  ns  at  Zs  is 

J^l 

E(n  at  Z/np  at  Zp^  ns  at  Zs)  =  ~]nF(n  at  Z/np  at  Zp^  n2  at  Z2) 

n=ns 


=  ng  +  (ni  -  ng) 


R(Z)  -  R(Z2) 

RrZi)  -  R(Z2) 


(8) 


where  N>n2_>n>n2>0  and  0  <  Z2_  <  Z  <  Z2 .  Equation  (8)  shows  that  the 
average  number  of  fibers  surviving  in  a  bundle  passing  through  n^^  at  Z-j^  and 
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rig  at  Zg  is  linear  in  R(Z).  Therefore,  the 
average  path  appears  as  a  straight  line  connecting 
the  points  ni  at  Zi  and  ng  at  Zg  on  a 
plot  of  n  against  R( Z ) . 


CLASSIFICATION  OF  PATHS 


In  the  calculations  in  this  section,  the 


n,  the 
average 


for  paths 
mum  value 


paths  will  be  grouped  in  accordance  with 
number  of  fibers  surviving  at  Z,  and  the 
loads  for  each  class  of  paths  shall  be  found  as  a 
function  of  Z  or  R( Z ) .  It  shall  be  shown  that 
the  average  load  for  jaths  1  (fig.  2)  having  more 
fibers  surviving  at  Z^  than  expected  reaches  a 
maximum  value  at  1=1,  whereag  the  average  load 
2  having  fewer  fibers  surviving  than  expected  at  Z  reaches  a  maxi¬ 


load  !,  R(0 
Figure  2.  -  Average  paths  of  bundles. 


at  Z  <  Z .  This  subdivision  of  paths  in  accordance  with^the  value  of 
n  carries  through  at  all  values  of  Z  so  that  the  average  load  L  at  any  Z 
for  bundles  having  a  given  n  may  be  determined.  The  subdivision  can  also  be 
regarded  as  a  first,  although  primitive,  step  in  the  approach  to  the  ideal 
wherein  the  probability  and  maximum  supported  load  are  determined  for  each 
possible  path,  so  that  an_averaging  of  the  maximum  loads  over  the  infinite  pop¬ 
ulation  of  paths  yields  S,  the  average  breaking  strength  of  the  bundle. 


In  accordance  with  eq_uation  ^8),  the  expectation  value  of  n  at  Z 
a  bundle  passing  through  n  at  Z  is 


for 


E(n  at  Z/n  at  Z)  =  < 


N  -  (N  -  n)  5^ 
B(Z) 

^  R(Z) 

.  R(0 


n  >  n,  Z  <  Z 

n  <  n,  Z  >  Z 


(9a) 

(9b) 


Figure  2  represents  the  situation  that  prevails.  Every  path  must  go  through 
the  initial  point  n  =  N,  R(z)  =1  and  the  terminal  point  n  =  0,  R(Z)  =  0. 

If  no  other  point  along  the  path  is  specified,  the  path  averaged  over  the 
bundles  is  given  by  curve  0  of  figure  2,  that  is,  by  the  straight  line  connecting 
the  initial  and  terminal  points.  If  a  point  is  specified  along  the  path  between 
the  initial  and  terminal  points,  two-segmented  lines  result  that  represent  the 
average  paths  of  bundles  going  through  the  additional  point  specified  (curves  1 
and  2  of  fig.  2).  The  additional  point  specified  for^curve  1  is  n^  at  Z, 
whereas  the  additional  point  for  curve  2  is  ng  at  Z . 

When  the  load  per  fiber  is  Z^  the  average  total  load  supported  by  a 
bundle  which  goes  through  n  at  Z  is 

L(z/n  at  Z)  =  ZE(n  at  Z/n  at  Z)  (lO) 
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MAlXIMM  loads  on  paths 


^The  maximum  load  that  occurs  on  an  average  path  when  a  particular  n 
at  I  is  specified  can  be  determined  as  follows.  Let  G( I )  =  ZR(Z);  then 
equations  (l)  and  (3)  can  be  written 


L(Z)  =  NG(Z) 


(la) 


G'(Z)  =  0 


(3a) 


The  assumption  is  also  made  that  the  form  of  R( ^ )  is  such  that  L( Z )  exhibits 
one  and  only  one  peak^  that  peak  occurring  at  Z,  JVIathematically  this  is 
equivalent  to  requiring  that 


G'(Z)  >  0 

z  <  z 

G'(Z)  <  0 

A 

(Ha) 

(lib) 


Horn  equations  (9)  and  (lO)^ 


dL( Z/n  at  Z ) 
dZ 


n  -  NR(Z)  +  (N  -  n)G'(Z) 
B(Z) 


[3(1) 


I  <  I 

I  >  t 


(12a) 


(12b) 


Equation  (l2a)  shows  that,  for  values  of  n  >  KR(Z),  L(z/n  at  Z)  is 
always  increasing  up  to  1  =  1.  At  Z  =  Z,  the  slope  changes  sign  (eq.  (l2b)) 
and  L(z/n  at  Z)  decreases  as  Z  in£reases  when  Z  >  Z.  Hence,  when 
n  >  NR(z),  the  maximum  value  of  L(z/n  at  Z)  occurs  at  Z  =  Z  or 


L(z/n  at  Z)^g^x  =  nZ  im(Z)  <  n  <  N 


(13) 


Again,  by  equations  ^12),  for  values  of  n  <  NE(z),  L(z/n  at  Z )  reaches 
a  maximum  when  ^  <  Z,  where  Ag  satisfies  the  equation 


(14) 


With  6£  denoting  the  right  side  of  equation  (14)^  the  following  quantities 
can  be  expanded  in  powers  of  5g,  where  the  coefficients  of  the  powers  of 
are  functions  of  the  derivatives  of  G( Z )  evaluated  at  Z : 


^5  p,2 

2G| 


G(A£)  =  G(z)  +  ^  4 


1  ^ 

3  r.3 
*^2 


(15) 

(16) 


7 


(17) 


L( ye  at  ?)  =  nl  -  i.  ai .  I  .  ^)6i  t  0(4) 

2 

In  equations  (15)  to  (17)  the  symbol  G-q  represents  the  n^^  derivative  of 
G{l)  with  respect  to  I  evaluated  at  1=1.  Note  that  G2  <  0. 

Equations  (l3)  and  (l?)  show  that  to  first  order  in  6£,  when  the  paths 
are  classified  in  accordance  with  n  at  I,  Daniels'  result  holds,  that  is, 
the  breaking  strength  of  the  bundle  is  only  dependent  on  the  number  of  sui*vivors 
at  Z . 

Equation  (l7)  may  be  employed  to  find  a  larger  lower  bound  for  S,  the 
average  breaking  strength,  than  that  given  by  equation  (2): 

N 

S  >  P(n  at  Z)L(A-/n  at  Z)  >  NZR(Z) 


[M(Z)] 

I  2 


If  the  variable  is  changed  from  n  to  y,  where  y  is  of  standard  measure 
and  is  defined  by 


-■  KR(Z)-n 


the  frequency  function  f(y)  can  be  expanded  in  terms  of  Hermite  polynomials 
Hn(y)  hy  the  method  of  Kendall  and  Stuart  (ref.  3): 


4=  e-y^/2 


fl  +  I  ^  (y3  -  5y)  +  .  .  .1 
L  Vhe(i)b(z)  J 


Also  6^  can  be  written  as  a  series  in  powers  of  y: 


(-1)^ 


^ilT 

Liro(z)J 
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Substituting  in  equation  (l8)  yields 


S  >  IltR(t)  + 


00 

-  f  — 


:-y^/2 


[l  +  i  1' )  „ ,  (y3  -  3y)  +  .  .  .1 

L  ®  V^lffi(Z)B(Z)  J 


[■ 


i  R(Z)  2  +  +  l\rR(Z)p/^  _zi 

G2  bTzT  \3G3  G2;[b(z)J  ^ 


•  •  • 


>  NZR(Z) 


i_  Tl  +  2.  2R(^)  -  1 


^2  B(Z)  L  3  ■y^2itM(Z)B(Z) 


■] 


m 


^5  ,  1 

3G|  ^2 


h=r-  NPT 
BitN  Lr(Z)J 


./Z 


■^itN 


.]+  0(N“^)  (22) 


The  variance  of  the  breaking  load  based  on  the  previous  approximations  is  given 
as 


0^ 


M(Z)B(Z)(Z)2  + 


2R(Z)  -  1 

V'kr(z)b(z) 


+  • 


+  0(1) 

(23) 


CONCLUDING  REMARKS 

A  method  of  determining  a  lover  bound  for  the  average  breaking  strength  of 
a  bundle  of  classical  fibers  has  been  presented  vhich  involves  the  subdivision 
of  the  bundles  in  accordance  vith  the  number  of  unbroken  fibers  at  a  given 
value  of  load  per  fiber.  The  method  also  yields  the  approximate  statistical 
distribution  of  the  breaking  strengths  so  that  the  variance  of  strengths  may  be 
calculated.  The  result  agrees  vith  the  literature. 


Levis  Research  Center^ 

National  Aeronautics  and  Space  Administration, 
Cleveland,  Ohio,  September  13,  1965. 
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aeronautical  and  space  activities  of  the  United  States  shall  he 
conducted  so  as  to  contribute  ,,,  to  the  expansion  of  human  knowl¬ 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  disseminatio72 
of  information  concerning  its  activities  and  the  results  thereof.” 

— ^National  Aeronautics  and  Space  Act  of  1958 
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